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Abstract 

Recently it was shown, within the mass insertion approximation (MIA), 
that the gluino mediated flavor changing neutral current (FCNC) interactions 
can saturate both ex and Re(e'/eftr) even for the real Cabibbo-Kobayashi- 
Maskawa (CKM) matrix through a single CP violating parameter (#f 2 )iL- In 
this work, we extend our previous analysis to the nonvanishing KM phase, 
and to the effective SUSY model where the MIA is no longer a good ap- 
proximation. In our model, €k and B° — B° mixing can receive significant 
SUSY contributions. Therefore the usual constraints on p and n in the 
Wolfenstein parametrization of the CKM matrix can be relaxed except for 
the iKtfe/^cbl from the semileptonic b —* u transition. This affects K — > ttvv 
and Kl — ► p + p~ indirectly, even if they are not directly induced by gluino- 
mediated FCNC, while Kl — > 7r°e + e _ is influenced both directly and indi- 
rectly. Differences between our model and other SUSY models enhancing 
Re(e'/eR-) are discussed in detail. 
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I. INTRODUCTION 



For many years, CP violating phenomena was observed only in Kl — > 2-7T M, which could 
be attributed to K° -K® mixing [AS = 2). The CP violating parameter in the K° — K° 
mixing, ex, has been accurately measured : ex = e i7r//4 (2.271 ± 0.017) x 10~ 3 f2|. However, 
the recent observation of Re(e'/e^) by KTeV collaboration, Re(e'/e^) = (28.0 ±4.1) x 10~ 4 
0, nicely confirms the earlier NA31 experiment [|J Re(e'/e^) = (23 ± 7) x 10~ 4 . Including 
another new data from NA48 collaboration, (14.0±4.3) x 10~ 4 [|5|, the current world average 
for Re(e'/e^) is 

Re(e'/e K ) = (19.2 ± 2.4) x 10" 4 . (1) 

This nonvanishing number indicates unambiguously the existence of CP violation in the 
decay amplitude (AS* = 1), and the original form of superweak model is excluded. Along 
with the measurements of sin 2/3 from time- dependent asymmetry in B° — > J/ipKs at B 
factories ||, we are now in the blooming era for detailed studies of CP violations in K and 
B systems. However, in this article, we will consider CP violations in the K system only for 
the reason to be explained later. 

Two parameters ex and Re(e'/e^) that quantify CP violations in the kaon system can be 
accommodated by the Kobayashi-Maskawa (KM) phase in the Glashow-Salam-Weinberg's 
standard model (SM). However, the SM predictions for the latter is rather uncertain because 
of large uncertainties in the nonperturbative matrix elements and the strange quark mass. 
There are significant variations in theoretical predictions [[/|-|TT|j, which indicate the current 



status of our understanding of Re(e'/e^). In this work, we will follow closely the work by 
Bosch et al. [0] in calculating Re(e' /ex) and other rare kaon decays so that we can compare 
the numerical results with other results more easily. 

Although the SM can accommodate a bulk of the observed Ke(e'/ex) within certain 
theoretical uncertainties, it is also interesting to speculate possible new physics contributions 
to CP violations in the kaon system. Among various scenarios beyond the SM, the minimal 
supersymmetric standard model (MSSM) and its various extensions are most promising 
candidates. Thus it is quite natural to consider generic SUSY contributions to ex and 
Re(e'/ex)- In this regard, two of us (SB and PK) showed that the flavor preserving A t and 
/i phases in the effective SUSY models alone cannot generate large enough ex |T^-4T4| or 



new phase shift in B° — B° mixing beyond that coming from the KM phase |15| . For e /ex 



there have been a lot of works done after the KTeV announcement fll6l-p9H . Those works in 
the supersymmetric models can be divided into three categories : 



the enhanced sdg vertex p9|-p7[] 



the enhanced sdZ vertex p6| , |29" l 



the enhanced EW penguins |E8 



Models in the first category are divided further into two parts, depending on the origin of 
the flavor and CP violations : either from the flavor changing squark masses in the LL sector 
or the flavor changing LR sector due to nontrivial flavor structures in the trilinear A terms 
which may arise naturally in string inspired models with D branes. Our model presented 
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in Ref. [23] is unique compared to other models with enhanced sdg vertex, since SUSY 



CP violations in our model affects not only the Re(e'/ex) but also the e# by significant 
amounts, contrary to other models in this category. In fact, it was shown that both tx 
and Re(e'/ejc) can be saturated by a single complex number (5f 2 )LLi if |/^tan/?| is large 
enough ~ O(10) TeV so that the induced (Si 2 )lr through the double mass insertion can be 
O(10~ 5 ). Furthermore, there is no conflict with the neutron electric dipole moment (EDM) 
in our model. There is another work showing that all the CP violations observed so far 
can be accommodated by a single CP violating phase in the gluino mass parameter M- g 
They assumed a very specific flavor structure in the trilinear A couplings in order to 



generate large e . On the other hand, we assume a specific flavor structure in the left squark 
mass matrix only, and the flavor structure and CP violating phase of a trilinear A coupling 
is totally irrelevant in our case. 

In generic SUSY models, there can be potentially large contributions to flavor changing 
neutral currents (FCNC) and CP violating processes from gluino-squark mediations, which is 
termed as SUSY flavor/CP problems [pOj] . Since these effects are due to strong interactions, 



it may be parametrically large compared to the charged Higgs and chargino contributions, 
unless the squarks and gluinos are too heavy. There are basically three ways to solve SUSY 
flavor problems : (i) universality of squark masses, (ii) alignment of squark and quark mass 
matrices due to some flavor symmetries, and (iii) decoupling of the lst/2nd generation 
squarks. 

In order to study the gluino mediated flavor changing phenomena within the frameworks 
of approximate universality or alignment, it is convenient to use the so-called mass insertion 
approximation (MIA) |3]J . In this approximation, one works in the super CKM basis where 
the quark mass matrices are diagonal. The quark-squark- gluino vertex is flavor diagonal 
in the MIA, and the flavor /chirality mixing occur through the insertion of (Sf^AB, where 
i, j = 1, 2, 3 and A, B = L, R denote the flavors of the squarks under consideration and the 
chiralities of their superpartners. The superscript denotes that the down type squark mass 
matrix is involved. The parameters (<5^)ab characterize the size of the gluino- mediated flavor 
changing amplitudes, and they may be CP violating complex numbers, in general. In the 
following, diagrams involving charged Higgs, chargino and neutralino will be ignored, since 
they are suppressed by a w /a s compare to the gluino-squark loops unless gluino/squarks 
are very heavy. This should be a good starting point for studying the SUSY FCNC/CP 
problems. 

In addition to the MIA, we also consider the vertex mixing (VM) model which is a good 
approximation if the third generation sfermions are lighter than the lst/2nd generation 
sfermions as in the effective SUSY models E|. In this case, the SUSY FCNC and CP 
problems are solved by the decoupling (and some degeneracy in the lst/2nd generation 
sfermion masses), and we keep only the left/right sbottoms and gluino contributions to the 



AS = 1, 2 effective Hamiltonians. The related works in B physics can be found in Refs. [33 



Although SUSY CP violations can saturate both ex and Re(e'/ex), it would be unnatural 
to assume that 5km vanishes (or equals 7r), since there is no symmetry principle supporting 
such assumption []. In this work, we present more detailed analysis of our model within both 



We can impose CP symmetry to be broken only softly. This case is already included in our model 
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MIA and VM cases with nonvanishing 5km- Since the ex and B° — B° mixing are affected 
by SUSY contributions significantly in our model, the usual constraints on (p, rj) from these 
two quantities need not be held any more. Only the constraint from b —>■ u semileptonic 
decay would be valid. Therefore, when we vary the (p, rj) , we will impose this constraint 
only, and the apex of the unitarity triangle (UT) can be anywhere in the (p, 77) plane, not 
only in the first quadrangle. This is in sharp contrast to other analysis within SUSY models 
where new physics contributes significantly only to Re(e'/e#) or the analysis based on the 
MSSM with minimal flavor violation. There is an indication that the unitarity triangle does 
not collapse (namely a, /3, 7 do not vanish) just from determining the triangle only from the 
unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [35| : 



a = 19° - 142°, /3 = 6°- 31°, 7 = 28° - 152°. (2) 

Also there are some indications that charmless nonleptonic 2-body decays of B mesons seem 
to prefer cos 7 < [56,37], although theoretical uncertainties involved are of vastly different 
degrees. Considering all these points, we consider all the possibility, although 7 ~ tt may 
be unrealistic. 

Assuming that the sk and B° — B° mixing cannot be used to constrain the CKM elements, 
we study several rare K decays such as K —>■ ttuu, K — > irl + l~ and K — ► p + p~ including 
the squark-gluino contributions, and compared with other recent works on Re(e /ex)- We 
do not consider the CP asymmetries in K — > 3tt, K — > 7r7r7 and hyperon decays. These 
observables were discussed in the literatures within the context of the enhanced sdg vertex as 
the origin of large e /ek- In our model, we have another source of FCNC and CP violation, 
namely (5( 2 )ll which induces AS = 1 four-quark operators. The detailed study of these 
effects to K — > 3tt, K — > mry and hyperon decays are beyond the scope of the present work. 
Also these processes will suffer from large hadronic uncertainties as other nonleptonic kaon 
decays in chiral perturbation theory, unlike K — ► tti/P. 

The contents of this paper are organized as follows. In the Sec. II A and B, we review 
briefly the AS = 2 and AS = 1 effective Hamiltonians and their relations to €k and e /ex- 
In Sec. II C and D, we recapitulate the basic formulae for neutron EDM, and branching 
ratios and CP asymmetries for K — > tti/P, K — ► irl + l~ and K — > p + p~ . The numerical 
inputs for various parameters in our model are summarized in Sec. II E. In Sec. Ill and IV, 
we discuss the CP violation phenomenology in the kaon sector for the MIA and the VM 
cases, respectively. Finally we summarize in Sec. V. 

II. EFFECTIVE HAMILTONIANS FOR KAON PHYSICS 

In order to discuss ex, e'/eg- and other rare kaon decays within the SM and SUSY models, 
we first construct the effective Hamiltonians for AS = 2 and AS = 1. We obtain the Wilson 
coefficients of the effective theory after matching it to the full theory at 0(my/). In this 
section, we define the operator basis for these effective Hamiltonians and the recipes for the 
renormalization group (RG) running from the weak scale to the lower energy scale p ~ m c . 



with 5km = or tt. See Refs. [34| for more discussions 
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The relevant Wilson coefficents for the AS = 1 and AS = 2 effective Hamiltonians are 
explicitly given in the Sees. Ill and IV, where we describe the MI and the VM approximations 
for the CP and flavor-violating effects of SUSY in detail. 



A. Effective Hamiltonian for AS* = 2 and ck 



The effective Hamiltonian describing the AS = 2 K° — K° mixing in the SM and the 
MSSM (with only the gluino contributions included) can be written as 



H, 



AS=2 
eff 



where 38.39 



TT 



AS=2 
SM 



G 2 p 



-m 



w 



\*M SM) s (x c ) + xfvr vi> s (xt) + 2\* c \;4* m> S (x c ,x t ) 



ttAS=2 i ttAS=2 
n SM "r -"SUSY 



(3) 



(SM) , 



x [« s (/x)]- 2/9 ^(1-75)^7^(1-75)/ + h.c, 



AS=2 
SUSY 



I 



216 m 2 



( 2 )/ 



(4) 
(5) 



\i=l 



i=l 



with operators defined to be 



Q 1 = a 2 s (fi) J q 7m(1 - 7s)s a JV(l - 75)/, 
Q 2 = a 2 s (fi) d a (l - 75)^^(1 - 75)/, 
Q 3 = a 2 s {fx) d a (l - 75)^(1 - l,)s a , 
Q 4 = a 2 s {fi) d a (l - 75)^^(1 + 75)/, 
Q 5 = a 2 (fi) d a (l - 75)^(1 + l,)s a . 



(6) 



Here, m is the average squark mass in MIA, or the sbottom mass in VM. a and (3 are color 
indices. The operators Qi are obtained from Qi by interchanging L and R. The explicit 
form of the loop function, So (2, 2/) ; can De found in Ref. |38 . 



We have separated the effective Hamiltonian into two parts : one coming from the SM 
and the other from SUSY sector, because the former is of 0(a w ), whereas the latter is 
of 0(a s ). We separately do RG running of SUSY and SM contributions to get the Wilson 
coefficients for the kaon physics at 0(m c ). In the SM part, we use the well-known anomalous 
dimesion matrices for RG running which can be found in Ref. [38|. For the RG running of 



SUSY part, we use a modified method advocated in in which new effective operators 
are defined including the a" factor coming from the gluino-squark loop diagrams. 

After defining the new operator set, we can apply the RG procedure of SM to gluino 
mediated contributions: matching calculation for Wilson coefficents to order a° and RG 
evolution by using the anomalous dimesion matrix to order a] for leading log approximations 
(LLA). The calculation of matrix elements is obviously modified by the additional factors 
of new operators. 

The anomalous dimension matrix in our basis at leading order is given by 



7 



(0) 



7<°> +4(3 ( J > 1 



(/) 



(7) 



5 



where (3$ 



'11 — 2//3) with the effective flavor number of / and 7^°- ) can be obtained in 
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(o) 



/4 0\ 

-28/3 4/3 

16/3 23/3 

-16 

V -6 J 



(8) 



The Wilson coefficients at 0(m c ) running down from those at 0(mw) read 

CfVc 



with 



and 



t(2). 



d 2) fm, 



r] 22 C ( 2 2) (m w ) + T^Cf^m^), 
7732 C| 2) (m w ) + ^33C3 2) (^w), 



V5 



a. 






(™c) 


a. 




a s 


(m c ) 


a, 




a { 


(™c) 


Q', 




a { 


(j™c) 


a s 


(m b ) 


a. 


(m c ) 



56/25 



1.419 



2.661 



26/25 



53/25 



a 



: (m w )\ 



52/23 



a s (m b ) J 



a 8 (m w 



1.369 



V22 = 0.983772 + 0.017//3, 
7/32 = -0.064r/ 2 + 0.064r/3, 



v a s (m b ) J 

I \ \ 2-718 

a s {m w ) \ 
a s {m b ) J 

/ a s (m w ) \ 22/23 
v a s (m fc ) / 

v a s (m fe ) j 



7/23 = -0.258//2 + 0.25 87/3, 

7/33 = 0.0177/ 2 + 0.983773- 



(9) 



(10) 



(11) 



In order to estimate AM^ and €k, it is necessary to know the hadronic matrix elements: 
(K°\Qi(fi)\K°) for fi ~ m c . In the present work, we use the following expressions [52, |3 
which are sufficient for our purpose in the LLA. 



(K°\QiM\K°) = -at l M m 2 K f K BM, 

(K°\QMm = -| o2(a*) f pfg n 

3 \m s (n) + m d {ji) i 



m 2 K fx B 2^), 



6 



3 \m s (n) + m d (fjL) 



{K°\Q^)\K°) = * 2 An) 



1 



1 + 



\ 2" 



\m s (/j,) + m d {ij) i 

2 / mt: 

3 lm s (/i) + m d (//) 



m2 Kf 2 K B i(v)i 
™k f 2 K B M- 



(12) 



The bag parameters 5's are given by 



B x (n) = 0.60, 
B 2 (» = 0.66, 
B 3 (n) = 1.05, 
B 4 (fi) = 1.03, 
S 6 (/i) = 0.73, 



(13) 



at fi — 2 GeV ( We use these values at /i = m c = 1.3GeV). 

The CP violation in the K° — K° mixing is characterized by a complex parameter ex, 
which can be obtained from the AS = 2 effective Hamiltonian by the following formula: 



r A = [ImM 12 + 2£ReM 12 ] 



(14) 



V2AM K 
where M i2 is defined as 

2m K M 12 = (K°\H* s = 2 \Kt). (15) 
Remember that £ = ImAo/KeAo is very small and thus is neglected in our calculation. 

B. Effective Hamiltonian for AS = 1 and Re (e ' /ex) 

The effective Hamiltonian describing the AS = 1 kaon decays receive contributions from 
the SM and reads (with only the gluino contributions included): 



r7-A5=l _ ttAS=1 , ttAS=1 
D cff — n SM -"SUSY J 



where 
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s ,, - c[o; + ^o 2 + E<^ + h.c. 



(16) 



(17) 



#s A usy = CiOi + C 2 2 + ^(CiOi + OA) + E E E CflxOiL + h - c - ( 18 ) 

i=3 A=± -B=7,g X=s,<? 

The explicit form of 0- can be refered in Ref. |38j and new operators Oj's are defined as 

O 4 (a0 = a s 2 (/i) o;(a*). (19) 



i=3 
10 



7 



Oi are obtained from O, by interchanging 1 — 75 and I+75. The 10 x 10 anomalous dimension 
matrix for RG evolution in these new four-quark operators is defined in a similar way to 
that of AS = 2 case given in Eq. (|7|). 

The last part of Eq.fllTf) is for the magnetic and chromomagnetic operators in which the 
subindex B denotes whether the chirality flips on the strange quark or the gluino line: 



O 



19 
± 



' [^^(1 + 75)^^^(1-75)5 



32tt 2 



ga s (fx) 



O 



± 

99 



327T 2 



da^t a G^(l + l5 )s ± da^t a G%(l - 75)* 
da^t a G« u (l + l5 )s ± da^t a G a ^(l - l5 )s 



The anomalous dimension matices for RG evolution of the above equations are given by 



(°) _ ^(o) 



,0 ' -= 7 (0) +2/?i /} 1, 



32/3 \ 
-32/9 28/3 ) ' 



(20) 
(21) 

(22) 



The Wilson coefficients of magnetic operators at 0(m c ) are 



Cg(g, s )(m c ) 



with 



a 

m 



and 



a s [m b ) 
a s {m c ) i 



<y s {m b ) 
a s {m c ) i 



29/25 



41/25 



' a s (m w ) \ 
a s (m b ) J 



' a s (m w ) \ 
v a s (m b ) J 



27/23 



V2 



, x \ 27/25 / 

aJm b )\ a s 



39/23 



V°2 



ac s [m c 



' a s (m b ) \ 
a s {m c ) ) 



a s (m b ) J 



25/23 



(23) 



(24) 



39/25 



QL s {mw) 
a s (m b ) 



37/23 



(25) 



The hadronic matrix element of a four quark operator in Eq. (|19|) is obtained by multi- 
plying the corresponding one of SM in Ref. |38| by a^(m c ). For the magnetic operators, the 
matrix elements at 0{m c ) can be expressed as pB] 



(O gs )o = m s (O aS ) = m s a s {{i) 



99' 



3 11 M^Ba, <£*>„ = 0, 



. 2167T 2 F3 

where the subscript means AI = components of the matrix elements. 



(26) 
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The parameter Re(e'/e_Rr), which characterizes the direct CP violation in the decay am- 
plitude of K — > 7T7T, is expressed as following in terms of the AS = 1 effective Hamiltonian: 

Hi) ^7^? Im(q) K (1 >') (27) 

where uj = Re^/ReAo = 0.045, and fl v + v ' = 0.25 ± 0.05 represents the isospin breaking 
effects. 



C. Effective Hamiltonian for the neutron EDM 



The electric dipole moment of a fermion / can be described by the effective Hamiltonian 



4l| 



H™ = ]T Cf m Of m , (28) 



i=l 



and its basis operators 



Of m = —fa^-rsfF^ 

Of™ = —fatcG^G^G^^. (29) 

With this Hamiltonian, we calculate the (chromo)EDM of quarks and gluons and compare 
the result to the experimental value to constrain possible parameter space. The current 
experimental upper bound on the neutron EDM is [55] 

|d n | <6.3 x 10" 26 ecm. (30) 



D. Rare kaon decays : K — » iwv, ttI + 1 , n + n 

SUSY contributions to various rare kaon decays through the (chromo) magnetic and en- 
hanced sdZ operators have been discussed comprehensively by Buras et al. |26|], and we 



heavily lean on formulae for branching ratios for rare kaon decays given in their paper. In 
the numerical analysis, however, there are two major different points to be emphasized : 

• In our case, we ignore the enhanced sdZ operator, which may be important only in 
a restricted parameter space, as already noticed in Ref. [F2IJ. Therefore, our model is 
precisely the Scenario A therein. 

• When we vary CKM matrix elements X t , we impose constraints coming from semilep- 
tonic decays only, since these are fairly insensitive to new physics beyond the SM, 
especially on the SUSY contributions we are considering here. On the other hand, 
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€k and -Brf( s ) — ^d(s) m i xm g ma y be affected by new physics. Indeed gluino mediated 
FCNC present in our model can affect both of these observables in a significant man- 
ner, unlike those considered in most recent works related with e /ex in the framework 
of SUSY. Since we consider the kaon sector only and the B meson sector is indepen- 
dent of the kaon sector in the MIA, we do not impose the B° — B° mixing in this work. 
Therefore, the weak phase 7 = 03 can be anywhere between and 2ir. 

Also {Sf 2 )LL would contribute to s — > dqq penguin operators, thereby modifying nonlep- 
tonic kaon decays (and CP violations therein) such as K — > 3tt and K — > irrcy and hyperon 
decays. The effects of (5f 2 )LR 011 these decays have been considered in the literatures, but 
not those of (Sf 2 )LL- Since these will involve some hadronic uncertainties compared to rare 
kaon decays we consider here (such as K —>■ ttvv), we will not include these nonleptonic 
kaon and hyperon decays in this work. 

In the SM, K — > txvv is affected by the KM phase. In fact the neutral mode Ki — > tx°vv 
is purely CP violating so that its branching ratio would vanish in the limit of exact CP 
symmetry. In our model with two sources of CP violations, only 5km affects these decays, 
since s —>■ dg vertex does not contribute to them. 



Using the formulae given in Ref. [26| we have the results 



B(K L -> 7T°UU) 
B(K + -> 7T+Z/P) 



B° SM = 6.78 x 10- 4 (A ImA t ) 



B, 



SM 



1.55 x 10 



(AoImAi) 2 + (X ReA t + A c 



(31) 
(32) 



where X = Cq — AB Q = 1.52 is a combination of loop functions evaluated at rh t (m t ) = 
166 GeV, and A c = — (2.11 ± 0.30) x 10 -4 is the internal charm contribution. 

Unlike K — > irvv, the direct CP violating part of Kl — > ir°e + e~ is affected in our model, 
because the Wilson coefficients of 7 (§, s ) an d O7V = (^• S )(y-A)(ee)y are modified due to the 
new CP violating phase from the squark-gluino diagrams. As for other SUSY contributions, 
we work with the 'SUSY version' of operator 



n SUSY/ 
U 7V 



a a s (n) (ds)(y-A)(ee) v , 



(33) 



rather than O7V, i n order to do matching and running formally consistent with LLA. Under 
the definition, 



ImA+Z/ 7 



a s B T 



V2G 



F m K 



m s ImC+ + ImC+- 



(34) 



the branching ratio is given by 



B(K L 



7i°e + e )dir 



6.3 x 10" 



(lm\ t y 7V + ImA+y 7 + ^2na s {m c )\mC^ Y ^ (35) 



+ (lm\ t y 7A ) 



where, for 



325 MeV and m t = 170 GeV, IBS 
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yrv(n = 1-3 GeV) = 0.537 x 2tt, 
y 7A (fx = 1.3 GeV) = -0.700 x 2tt. 



The short distance part of the branching ratio of Kl — > H + fJ> within our model is the 
same as that in the SM [ 25 ]: 



B{K L -> ^ + /x-) S d = = 6.32 x 10- 3 (Y ReA 4 + A c ) 2 , (36) 
where A c = —(6.54 ± 0.60) x 10 -5 is the charm contribution. 

E. Input parameters 

Since there may be large SUSY contributions to e# as well as to B° — B° in our model, 
one cannot use the constraints on the CKM matrix elements coming from these observables. 
We impose only the constraint coming from the semilpetonic decays of mesons. Therefore 
we have [[J 

\V ub /V cb \ =0.08 ±0.02 (37) 

This leaves the phase of V ub undetermined so that the SM contributions to various observ- 
ables are not fully determined. In the following numerical analysis, we sample data fixing 
the phase of V ub to some selected values of 0°, 60°, 120°... so on up to 360°. 
We list numerical input values used throughout our analysis in Table |. 



III. MASS INSERTION APPROXIMATION CASE 

A. Approximate flavor symmetry and (Sf 2 )AB 

One way to solve the gluino-mediated FCNC problems in general SUSY models is to 
assume that the squark masses are approximately diagonal and degenerate in the basis 
where quark mass matrices are diagonal (the so-called super CKM basis). 

{rh%)AB = {m 2 )AB^ij + (Sm 2 i:j ) AB . 

In this case, one can expand the squark propagator around the common squark mass matrix 
{fh 2 )AB^iji and the rest can be cast into interaction Hamiltonians. As long as 5m 2 is small 
compared to the common squark mass m 2 , one can keep the lowest order perturbations 
in 5m 2 relevant to the processes we are interested in. In many cases, just a single mass 
insertion along a single squark propagator is sufficient, but double mass insertion along a 
single squark propagator should be included in some cases. In fact, it is the case for e /ek 
in the relatively large /xtan/3 region f23fl . In this case, smallness of s quark mass can be 
compensated by a large //tan/3 (which has been ignored in most previous literatures) so 
that an appreciable amount of Sr — §l may be possible. 

The typical size of (5£)ab = (8 m2 j) ab / ™ 2 should be fairly small in order to be consistent 
with AMk and ex- Theoretical understanding of such a small number constitutes the so- 
called SUSY flavor problems. This problem could be solved by approximate universality or 
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approximate flavor symmetry where both quarks and squarks are almost aligned in flavor 
space. 

The gluino-squark contributions to the AS = 1,2 effective Hamiltonian at the heavy 
SUSY particle scale have been calculated by several groups, and are encoded in the Wilson 
coefficients listed below. 



B. Wilson coefficient for AS = 2 

The Wilson coefficients for the AS = 2 effective Hamiltonian at the SUSY particle mass 
scales coming from the gluino contributions are given as follows: 

Cf = - 1 - (24x/ 6 (x) + 66/ 6 (a;)) (8 d 12 )\ L , 

r (2) _ r (2) _ r (2) _ r {2) _ p;(2) _ n 

cP-^xmiltfJlk)^. (38) 

The expressions for Cf } and cf ) are new, but their contributions to e K and Am K are 
negligible. The loop functions are given by the following expressions : 

6(1 + 3x) lux + x 3 - 9x 2 - 9x + 17 



f 6 (x) 



6(x- l) 5 



~. 6x(l + x) In a; - a; 3 - 9x 2 + 9x + 1 

/«(*) = » ( 39 ) 

197 + 25a; - 300a; 2 + 100a; 3 - 25a; 4 + 3a; 5 + 60(1 + 5a;) In a; 

where x = {rrig/ifi) 2 . The flavor conserving but chirality changing mass insertion parameter 
(^22) lr is defined as 

777 

0&)m = =?K - (40) 

777 

We neglect the terms proportional to (Sf 2 ) 2 LR, (Si 2 )rlj an d (^12) /ir- 



C. Wilson coefficient for AS = 1 

The Wilson coefficients for the AS = 1 effective Hamiltonian at the SUSY particle mass 
scales are given as follows : 

d = (for i = 1,2,7-10) 

°> = ISP - §*<*> - Ts Pl{x > - \ p ^) 



12 



C, = ^ (-Wx) + \b 2 {x) + l Pl (x) + |p 2 (x)) (St 2 ) 



' LL 



4m 2 V 3 v ' 3 v y 6 w 2 v V V "111 
Cj = Cj with replacement of LL by RR, % — 3, 4, 5, 6 (41) 



where 



Pi(x) 
P 2 (x) = x 
Pi(x) 



1 + 4a; — 5x 2 + 4xlogx + 2x 2 logx 



8(1 -x) 4 
5 — 4x — x 2 + 21ogx + 4xlogx 
2(1 -x) 4 

1 - 6x + 18x 2 - 10x 3 - 3x 4 + 12x 3 logx 
18(x- l) 5 



. * 7 - 18x + 9x 2 + 2x 3 + 31ogx - 9x 2 logx . . 

P ^ X) = 9(x^Tr (42) 

Bi and Pj are box and penguin diagram contributions respectively. 

(«) 

C t> = ~% {<) LL (Im,(x) - , (45) 

CS = § (™^) (§"><*> - 5 M »W) (46) 

where Qd = —1/3 and 



Mi (a) 



3-3x 2 + (l+4x + x 2 )lnx 



(x-1) 



5 



,,,, l + 9x-9x 2 -x 3 + 6x(l + x)lnx 

= 2(x^T)^ — (47) 

M 3 (x) = ^M 2 (x) 



M 4 (x) 



3 

1 + 9x + 9x 2 - 17x 3 + 6x 2 (3 + x) lnx 
12(x- l) 5 



Here again, the Wilson coefficients C^~ g and arising from double mass insertions are given 
for the first time. These two new terms distinguish our model from other models, when we 
discuss the processes Kl — > n°e + e~ and Kl — > fi + \i~ through s — > dj, and the process 
Re (e /ex) through s — > dg. The process Kl — > 7r°e + e~ is also affected by the s — > <i/ + /~ 
local interaction whose Wilson coefficient at the SUSY particle mass scale is 

CT Y ' = -^- 2 Pi(x)(5f 2 ) LL . (48) 
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D. EDM constraint 



The Wilson coefficients for the effective Hamiltonian for the neutron EDM in MIA are 
given by 

cfm = n [m £ Mx) _ 4 Bi(i) ^ f (50) 

where 

«)^ = il(^-/itan/?). (51) 

If we consider the loop functions to be order of one, and set both the squark and the 
gluino masses to 500 GeV, we get the limit on the flavor preserving mixing parameter, 

Im(^) Li? <O(10- 8 ), (52) 

from the constraint (|30|). The typical modulus of (5f x ) lr is O(10~ 4 ) for A^—fi tan (3 = lOTeV, 
and this limit constrains the phase of (S^lr below O(10~ 4 ). 

There is a correlation between phases of (<5fi) lr and (822) lr because they share a common 
term /xtan/3, and this correlation becomes stronger if we assume universality between 
and A s . This indirectly constrains the phase of {S^lr- But we can safely set its phase to 
zero while satisfying ex and Re(e' /ex) at the same time, as we will show in the numerical 
analysis. 



E. Numerical Results 

Now in Fig. [I], we show the plots of Re (e'/ e K ) as a function of the phase of (8f 2 )LL, where 
(^12) ll is parametrized as re 1 ^. We fixed A* — /itan/3 = 10 TeV. Different values and sign of 
A* — fi tan (3 yield similar results. For a given value of 0, r is determined in such a way that 
€k becomes the experimental value, and Re(e'/e^) is computed from those two parameters. 
As <fi varies from 180° to 360° for 7 = 0°, the SUSY contribution to Re(e'/ex) repeats its 
values in the region of <fi from 0° to 180°, with the sign flipped. Therefore it is clearly 
ruled out, and the results are not shown. But as the angle 7 varies, the SM contribution 
either increases or decreases Re^e'/ex), so that the invisible parts in 7 = graph float up in 
7 = 60° and 7 = 120° graphs, while some visible parts sink down in 7 = 240° and 7 = 300° 
graphs. Note that the functional relation between r and <fi governed by the ex constraint 
depends on the KM angle 7, and thus plots for nonzero 7 are not simple shifts of that for 
vanishing 7. 

Since there may be large contributions to K° — K° and B — B° mixing in our model, it 
is utterly important to measure three angles a,/3 and 7 of the unitarity triangle (UT) in a 
way indepedent of these neutral meson mixings, in addition to the more popular ways using 
B — > J/ipKs, 7TTT, so on. Also, it would be an interesting and important question how much 
portions of Re(e / /e^) come from the SM and new physics (SUSY gluino-squark loops in 



14 



this work). This question can be answered only if the angle 7 is determined by some other 
methods independent of B° — B° mixing, e.g., from B — > DK or B — > wK, etc.. Once the 
angle 7 is known, one can calculate the SM contribution, Re(e'/ex) SM - Then, the difference 
between the observed Re(e'/ex) and the SM would be the new physics contributions. 

In our model, K — > itvv decays are not directly affected by the gluino-squark loops. 
This is because neutrinos couple to neither of photons and gluinos, and thus the gluino- 
squark loop do not induce sdvv vertex. Still, the shape of the UT, namely the angle 7, 
can be different from the SM, because of potentially large contributions to (and also 
to B° — B° mixing). Therefore, gluino mediated FCNC will affect the branching ratios for 
K — > -kvv in an indirect way, even if the Wilson coefficients for s — > dvv remains the same 
as in the SM. In Fig. ^| (a) and (b), we show the Kl — > tt°vv branching ratio as a function of 
7 and its correlation with Re(e' /e^), respectively. There is no definite correlation between 
Kl — > tc°vv and Re(e'/ex) in our model, unlike the scenario in which the new physics affects 
Re(e'/ex) by a modified sdZ vertex. We note that one cannot make definite predictions for 
the branching ratio for this decay unless we know the angle 7 in the presence of new physics. 
The branching ratio may be vanishingly small or larger than the SM prediction by a factor 
of ;$ 2. Therefore, this decay once measured may provide us with invaluable informations 
on the value of the angle 7, if there is no significant new physics contribution to s -> dvv 
upto 4-fold ambiguities, which could be eliminated using information from other rare kaon 
decays and B meson decays. In Fig. |^ (c) and (d), we show the similar plots for K + — > tt + vv. 
Again, the measurement of the K + — > tt + vv branching ratio will help to determine 7. And 
the branching ratio for this decay may be larger than the SM predictions by a factor of ~ 2. 

In Figs. 3 (a) and (b), we show the branching ratio of B(Kl — > 7r°e + e~)dir as a function 
of 0, the argument of {Sf 2 )LL, and its nontrivial correlation with Re(e'/e^). This nontrivial 
correlation arises because electrons couple to a photon and also the s — > d'j vertex contributes 
both Kl — > 7T°e + e~ and Re(e'/ex)- In other words, the Wilson coefficient for s —>■ dj is 
modified by the gluino mediated FCNC, which is in turn closely related with the gluino 
mediated sdg operator. Also the Wilson coefficient of O-jy vertex is modified by the gluino 
mediated FCNC, which affects the direct CP violating part in Kl — > 7r°e + e~. This is in 
constrast with the K — > ttvv case which is not affected directly by gluino mediated FCNC. 
Therefore, the direct CP violating part of Kl — > 7r°e + e~ can be affected a lot, if there is a 
new CP violating phase in Sl — dL mixing both by a sing le MI of (5f 2 ) LL and a double MI of 
(^Ylr = (^12) ll x (#22) lr- Also this process is affected indirectly by SUSY gluino effects 
through their dependence on the KM angle 7. After all, the branching ratio of Kl — > ir°e + e~ 
can be vanishingly small or enhanced over the SM value by a factor of ~ 3. 

Finally, we consider the gluino mediated SUSY contributions to the short distance part 
of Kl — > Again this decay could be affected by gluino mediated FCNC both directly 

and indirectly : through the modified Wilson coefficient of sdj operator and modified values 
of 7 due to SUSY contributions to ex- However, unlike the direct CP violating part in 
Kl — > 7T°e + e~ , the modified sdj is irrelevant to Kl — > fi + {i~ because of current conservation 
associated with the muon current. Therefore it depends only on the angle 7 [Fig. 4 (a)], 
and its correlation with e'/e is rather trivial [Fig. 4 (b)]. Depending on the values of 7, the 
result can be reduced by a factor of ~ 3 or increased by a factor of ~ 2. 

Let us summarize the results on rare kaon decays within mass insertion approximation. 
We assumed there are significant gluino-squark loop contributions to ex and Re^'/e^) 
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through (Sf 2 )LL ~ O(10~ 2 — 1CT 3 ) and the {Sf 2 )LR ~ O(10~ 5 ) induced by a double mass 
insertion of flavor preserving Sr — sl followed by (Sf 2 )LL- Then these parameters directly 
affect s — > dg and s — > c/7 vertices, but not s — > dvv. Also the shape of the UT would 
be different from the SM case because SUSY loop contributes to e#. Therefore all the rare 
kaon decays (and also B decays as well, although we do not discuss this subject in any detail 
here) are affected indrectly through modified CKM elements, especially 7. For K — ► -nvv 
and the short distance part of Kl —>■ the branching ratios can be larger than the SM 

expectations by a factor of ~ 2 — 3, or can be smaller than the SM values. Especially the 
branching ratio for Kl — > 7i°uu can be vanishingly small if 7 = or tt. Their correlations 
with Re(e / /eft-) are rather trivial, because these decays are not affected directly by (Sf 2 )LL- 
On the other hand, the direct CP violating part of Kl — > vr°e + e _ is affected both directly 
and indirectly, and its correlation with Re(e'/ex) is nontrivial at all. One may expect much 
enhanced (or reduced) direct CP violations in K L — > 7r°e + e~ in our scenario, depending 
whether SUSY effects can interfere constructively (destructively) with the SM amplitude. 



IV. VERTEX MIXING CASE 

A. Flavor mixings in the quark-squark-gluino vertices 

Another way of solving the SUSY flavor and CP problem is to assume that the lst/2nd 
generation squarks are heavy and there are certain degrees of hierarchy in the qi — q~j — g 
mixing matrices, Wl and Wr (the decoupling scenario). Here Wl and Wr are analogous to 
the CKM matrix in the SM. Interactions among quarks, squarks, and gluinos in terms of 
their mass eigenstates are described by the Lagrangian: 2 

C = -V2g s (WL) ij q L i TRt a qi + ^3,^)^^ + ^, (53) 

where i and j are generation indices and color indices of (s) quarks have been omitted. Al- 
though squarks of each chirality are diagonalized into their (approximate) mass eigenstates, 
the remaining small mixing between squarks of different chiralities can be treated as a per- 
turbation except for the stop sector. 



B. Wilson coefficient for AS = 2 

The Wilson coefficients for AS = 2 effective Hamiltonian in the VM case are given by 
following expressions : 

Cf ) = +1 ( 24 / 4 (x) + 66/4^)) (Fn) 2 LL , 
Ci 2) = + 2 -^U(x)(F 12 )l L , 



2 Our definition of the mixing matrix W is different from that in Ref. [46 in that the signs of the 
two terms are opposite and both Wl and Wr are not complex-conjugated. 
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cf 
d 2 



c? 
c? 
c {2 



cf 



where 



~h{x) (F 12 ) 2 RL , 



+\ (504/ 4 (x) - 72/ 4 (x)) (F 12 ) LL (F 12 ) 
+± (24/ 4 (x) + 120/ 4 (x)) (F 12 ) LL (F 12 ) 
+i (24/ 4 (x) + 66/ 4 (x)) (F 12 )^, 



132 ~ 

-— h{x) (F 12 ) LR (F l2 ) 

RL ' 



4 
180 



/ 4 (x) (F 12 ) LR (F 12 ) 

RL > 



4 

+^/ 4 (x) (F 12 )^ 



f 4 (x) = X 



2x-2-(l + x) lux 



h{x) 



{x-lf 
x 2 — 1 — 2x Inx 



(x-1) 3 

(f 12 ) ab = (w A y 31 (w B ) 32 . 



(A,B — L or R) 



(54) 



(55) 



Here we assume that both left and right sbottom masses are equal to m, and define x 
{rrig/m) 2 . 



C. Wilson coefficient for AS = 1 

The Wilson coefficients for AS* = 1 effective Hamiltonian in the VM case are given by 
following expressions : 

d = (for i = 1,2,7 ~ 10), 
1 a 2 



C 3 
C 5 



C, 



± 

7s 



a 



■± 
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c. 



± 



99 



4 2m 2 
1 3a 2 



4 2m 2 
1 a 2 



(F 12 ) LL P(x), 
(F 12 ) LL P(x), 



4 2m 2 
1 3a 2 



(F 12 ) LL P(x), 



(F 12 ) LL P(x), 



4 2m 2 

Cj = Ci with replacement of LL by RR, i = 3, 4, 5, 6 
8nQ d 



3m 2 
87iQ d 



[(F 12 ) LL ±(F 12 ) RR ]P BE (x), 
rrib(Al — //tan/3) 



{F 12 ) 



LR- 



m 



3m 2 
2n 

- [(F 12 ) LL ±(F 12 ) RR ]P 2 (x), 



± (F 12 ) 



RL- 



m b (A b — ji* tan (3) 



m 



2ir 

2 



m 



m b (A* b - //tan/3) m b (A b - // tan/5) 

12) LR ~ ± {-t 1 12) RL 



m 



m 



(56) 

y/xP B i(x,x), (57) 
(58) 

y/ZP 3 (x,x), (59) 
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where the LR mixing parameter (_A£ — //tan/3) is real, and 

1 1 



Pf(x) 
Pb(x) 
Pfe{x) 
Pbe{x) 
Pfi(x) 
Pbi(x) 



Pbi(x.x) 



36 ( 


1 


-x) 4 


-1 




i 


36 


(1 


-x) 4 


1 




-1 


12 ( 


1 


-x) 4 


1 




1 



12(1 -x) 4 
1 —x 
2(1 -x) 3 
1 x 
2(l-x) 3 
1 



7x 3 - 36x 2 + 45x - 16 + (18x - 12) lnx 
llx 3 - 18x 2 + 9x - 2 - 6x 3 lnx 
x 3 — 6x 2 + 3x + 2 + 6x In x 
2x 3 + 3x 2 — 6x + 1 — 6x 2 In x 
Ax — x 2 — 3 — 2 In x 
1 — x 2 + 2x In x 



P(x) = Pp(x) - -Pb(x), 

3 1 
^2 (a?) = -Pfe(x) - -Pbe(x), 
2 6 

P 3 (x) = ~Pfi{x) - \pbi{x), 
I 6 

Pbi{x) - P B i(y) 



lim 



Ps(x, x) = lim 



x — y 
P 3 (x)-P 3 (y) 



'i— 1 x — y 



(60) 



We don't consider the box diagrams in vertex mixing, because it will be highly suppressed 
due to heavy lst/2nd generation squarks. 



The Wilson coefficient of Oyy SY/ is 



^SUSY/ 
°7V 



9m 2 



— P B (x)(F 12 ) LL . 



(61) 



D. EDM constraint 

The Wilson coefficients for the effective Hamiltonian for the neutron EDM in VM are 
given by 



^edm 2 ea s n rn- g m b Im [{F n ) LR {Aj - //tan/?)] , 
W = — — ^ F B i[x,x) 

3 7r m m z 

dm 1 g s a s m-g m b Im [(F n ) LR (A* b - //tan/?)] /3 1 x 

4 7T m 2 m 2 V2 6 



Basically, the constraint on lm.(6f 1 )LR in fl52|), applies to m fc Im [(F u ) LR (Al — //tan/?)] /m 2 
much the same way, and if we assume (A% — // tan (5) is real and is about 2 TeV, this limits 

lm(F n )LR^O(W- 7 ). (63) 
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Since we will set (Wr)3i to zero in the numerical analysis, there is no difficulty in satisfying 
this limit in our case, but it should be obeyed when other possible solutions are searched 
for. For example, if both (Wi) 3 i and (Wr) 3 i have absolute values of O(10~ 2 ), the typical 
order of magnitude of |(W / l)3i| in our analysis, then their phase difference should be less 
than £>(l(r 3 ). 

E. Numerical Results 

In MI case, the only key parameter responsible for CP violation was {Sf 2 )LL- Once 
its modulus is determined by ex constraint for a fixed 7, every CP violating observable 
could be computed as a function of its phase, Arg(5f 2 ) Li . Within the VM approach, on 
the contrary, the analysis becomes more complicated, because we have four independent 
complex parameters, (Wl)si, (Wl) 32 , {Wr)si, and (Wj^) 32 , that describe the flavor mixing 
in the q — q — g with various chiralities of quarks and their superpartners. Even if we impose 
the ex and the EDM constraints, there still remain 6 real parameters and it would be a 
formidable task to scan all of this 6 dimensional parameter space. In order to simplify the 
numerical analysis, we discarded some parameters and tried to find a particular solution 
rather than look for a completely general solution. That would not much distort the picture 
of VM approach and would suffice to illustrate a possible existence of appropriate parameter 
space in which gluino-squark loop contributions to €k and Re(e'/e^) are significant without 
conflict with electron/neutron EDM's. 

For definiteness, we set (1^)32 and (W 11)31 to zero, and all SUSY contributions are 
expressed in terms of a single complex number (-F 12 )lr = (Wl) 31 (Wr)^- The results for 
Re^e'/ex) for various values of the angle 7's are displayed in Figs. [|, as functions of 4>, where 
we use the following parameterization, (Fi2)lr = re 1 ^. As in the MIA case, we can have 
a large e'/e in certain region of with size ~ 0{1). In particular, all the CP violations in 
the kaon system could be accommodated in terms of CP violations from SUSY sector with 
vanishing KM angle, 7 = or ir, if Al — yutan/? = 0(2) TeV. Therefore, tan/3 could be 
considerably smaller than the MIA case because the flavor preserving LR mixing in the VM 
case is enhanced by m^jm s relative to that in the MIA. Also SUSY particle mass spectra 
will be quite different from those in the MIA case, although both cases can accommodate 
the CP violations in the kaon system. 

The decays K — > txvv and the short distance part of — > in the VM will 

be essentially the same as those in the MIA. There are no direct SUSY effects on these 
processes. SUSY effects are only indirect through the angle 7 from the SUSY contributions 
to €k- Therefore we do not show any plots for these processes. The 7 dependence are the 
same as the MIA [see Figs. 2 (a), (c) and Fg. 4 (a)]. 

On the contrary, the direct CP violating part of Kl — > 7r°e + e~ is affected directly by 
the SUSY loop contributions. In Fig. |6] (a), we show its branching ratio as a function of the 
angle for different values of 7. In Fig. ^ (b), we show its correlation with Re(e'/e^). It 
looks very similar to Fig. |^ (b) because = C~ and = C~ here as was in the MIA 
case. If we relax the assumption that (Wl) 32 = (Wr)3i = 0, however, then the correlation 
may differ from that in the MIA case in general because Re(e'/e^) mainly depends on C g 
while B(K L -> 7r°e + e-) on C+ 
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V. CONCLUSIONS 



In this paper, we discussed the gluino-squark loop contributions to CP violations in 
the kaon system both in the MIA and in the VM. Unlike many other models based on the 
enhanced sdg or sdZ vertex in the context of SUSY using nonuniversal flavor structures 
in the trilinear A couplings, our model modifies both the e and Re (e /e#), so that the 
phenomenological implications on other kaon decays substantially differ from other models. 
In particular the shape of unitarity triangle can be very different from the SM case, due to 
potentially large SUSY contributions to e K and B° — B° mixing. Only the constraint from 
the semileptonic b — > u transition is stable against a possible SUSY contributions. Therefore 
we varied the 7 from to 2n, then fixed the X t to get the right magnitude of Ek- In particular, 
one can accommodate both ex and Re(e /e#), even if the angle 7 is very close to or 7r, still 
consistent with the ranges obtained from the three generation unitairity relations, Eq. (2). 
In our model, the gluino-squark loop does modify s — > just like the sdg operator, but 
does not directly contribute to s — > dvv. Still the branching ratios for K — > tiuu can be 
different from the SM predictions by a large amount because of the modified 7. Similar 
effects occur in K L — > n e + e~ and K L — > fi + fi~ processes. If the shape of the triangle 
in MSSM is the same as in the SM case, then the deviations from the SM predictions are 
possible only in e /e^ and Kl — > 7t°e + e~, although the deviation in the latter case is very 
small. Also in our model, SUSY effects give additional contributions to AS = 1 penguin 
operators, so that we may anticipate interesting deviations of direct CP asymmetries in 
K — > 3tt and K — > tttt^ from the SM predictions. However, these observables are generically 
contaminated by hadronic uncertainties even in the chiral perturbation theory, and are left 
for the future study. 

One can also perform a similar analysis in the B meson sector including SUSY gluino- 
squark loop contributions to B° — B° mixing, and considering a specific flavor structure for 
the mass insertion parameters . Depending on the chiral structures of the mass insertion 
parameters, one could achieve sin 2(3j/^k s which is very different from the SM case. In this 
case, the results could be very different from the SM even if we assume that the UT in the 
MSSM is the same as in the SM. 

Our study indicates that the measurements of Kl — > tc^vv in the kaon sector is very 
important for constructing the unitarity triangle, especially if there are large new physics 
contributions to €k and B° — B° mixing. Also, in the B meson sector, it would be utterly 
important to construct the UT from the tree level processes which are immune to possible 
new physics contributions. For example, extracting 7 from B — > DK or other tree dominated 
decays will be very desirable ||48|| . The current global analysis on CKM matrix elements 



suggests that a unique UT emerges and the SM picture for flavor and CP violations in 
the quark sector seems to be a correct picture. But as we demonstrated with our specific 
model, this is no longer guaranteed in the presence of new physics which give significant 
contributions to ex and B° — B° mixing p9 [. 
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FIGURES 
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FIG. 1. The plots of e'/e versus the phase of (Sf 2 )LL in MI for six different values of 7 with 
A* — n tan /? = lOTeV. Graphs were drawn in the solid, the dashed, the dotted, and the dash-dotted 
lines for x = 0.3, 1.0, 2.0,4.0, respectively. 
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FIG. 2. The 7 dependences of B(K — > 7n/zv) ((a), (c)) and their correlations with e'/e in MI 
((b), (d)), for x = 1 with ^4* — /u tan/3 = 10 TeV. (b) The solid, the dashed, and the dotted lines 
correspond to 7 = 30°, 60°, 90°, respectively. The first two equally correspond to 7 = 150°, 120° as 
well. The branching ratio vanishes for 7 = 0° or 180°. (d) The solid, the dashed, the dotted, and 
the dash-dotted lines correspond to 7 = 0°,60°, 120°, 180°, respectively, and 360° minus them. 
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(a) (b) 

FIG. 3. The 4> dependence of B(Kl — > 7r°e + e~)di r ((a)) and its correlation with e'/e ((b)), 
in MI for x = 1 with ^4* — //tan/? = 10 TeV. (a) The solid, the dashed, the dotted, and the 
dash-dotted lines correspond to 7 = 0°, 180°; 60°; 240°, 300°; 120°; respectively, (b) The solid, the 
dashed, and the dotted lines correspond to 7 = 0°, 60°, 240°, respectively, and 180° minus them. 




(a) (b) 

FIG. 4. The 7 dependence of B(Kl — > /U + //~)sd (( a )) and its correlation with e'/e in MI ((b)), 
for x = 1 with A* — /itan/3 = 10 TeV. (b) The solid, the dashed, the dotted, and the dash-dotted 
lines correspond to 7 = 0°, 60°, 240°, 180°, respectively, and 360° minus them. 
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FIG. 5. The plots of e'/e versus the phase of (Fi2)lr in VM for six different values of 7 with 
Al — ft tan (3 = 2TeV. Graphs were drawn in the solid, the dashed, the dotted, and the dash-dotted 
lines for x = 0.5, 1.0, 2.0, 4.0, respectively. 



27 



-10 

10 

0.1 r- 

0.09 E- 

0.08 E- 

0.07 E- 

0.06 E- 

0.05 E- 

0.04 E- 

0.03 E- 

0.02 E- 

o.oi E- 





-10 



50 100 150 200 250 300 350 

<K°) 



10 

0.1 

0.09 
0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




e'/e (x 10" ) 



(a) (b) 

FIG. 6. The (j> dependence of B(Kl — > 7r e + e _ )di r ((a)) and its correlation with e'/e ((b)), in 
VM for x = 1 with yl* — \i tan (3 = 2 TeV. These curves are specific to the solution of parameter 
space drawn in Fig. || which restricts {Wl)z,2 and {Wr)^ to zero. The meanings of line patterns 
are the same as in Fig. ||. 
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TABLES 



m z 


91.2 GeV 


a s {m z ) 


0.118 


711 117 


8D 2 GeV 


f 


1 31 MpV 

1<J1 1*11/ V 


m t 


170 GeV 


Sk 


160 MeV 


m b 


4.4 GeV 


AM K 


3.51 x 10" 15 GeV 


m c 


1.3 GeV 


sin 2 9w 


0.23 


m d 


8 MeV 


\ £ k\ 


(2.266 ±0.023) x 10~ 3 


m n 


135 MeV 


(Re^o)cxp 


3.33 x 10~ 7 


m K 


498 MeV 


(ReA 2 )exp 


1.50 x 10~ 8 



TABLE I. Input values we used in the numerical analysis. 
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